Foreign exchange market fluctuations as random walk in demarcated complex plane 
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We show that time-dependent fluctuations {Aa;} in foreign exchange rates are accurately described 
by a random walk in a complex plane that is demarcated into the gain (+) and loss (— ) sectors. Ax 
is the outcome of A'^ random steps from the origin and |Aa:| is the square of the Euclidean distance 
of the final A'^-th step position. Sign of /S.x(t) is set by the A'^-th step location in the plane. The 
model explains not only the exponential statistics of the probability density of {Aa;} for G7 markets 
but also its observed asymmetry, and power-law dependent broadening with increasing time delay. 

PACS numbers: 89.65.Gh, 02.50.-r, 46.65.+g, 87.23. Ge 



Easier data access via the Internet and the widespread 
availabihty of powerful computers, have enabled many 
researchers not just in business and economics but also 
inphysics, applied mathematics, and engineering 0, |^ 
0, |a| , to investigate more deeply the dynamics of foreign 
exchange markets. Their efforts have led to new insights 
on the characteristics of forex rate fluctuations including 
the general behavior of their corresponding probability 
density function {pdf) 

Arguably, forex markets have a more immediate and 
direct impact on citizens than stock markets do. This 
is especially true in developing countries where only a 
minority holds stocks directly while a majority sells or 
purchases products and commodities on a daily basis. 
Economies that rely strongly on remittances of overseas 
contract workers (e.g. Philippines) or tourism (e.g. Thai- 
land), are also quite sensitive to forex rate instabilities. 
Hence, knowledge of the vagaries of the forex market, is 
important not only to policymakers and macroeconomic 
managers in government and banking sectors but also to 
individual businessmen and consumers. 

The exchange rate of a currency relative to another 
(usually the US dollar) is represented by a time-series 
of quotations {x{t)}, where time t = mAt, At = 
time delay, and index m — 1,2,...,M. Normally, the 
power spectrum of {x(t)} consists of one or two domi- 
nant low-frequency components and a continuous high- 
frequency band. The low-frequency components reveal 
the longterm behavior of {x{t)} and therefore the sound- 
ness of a national economy from the perspective of eco- 
nomic fundamentals. The high-frequency components 
are associated with fluctuations in {x(t)} that arise from 
complex interactions between market players over short 
time-scales (At < 2 days) [1,I3|- 

The statistics of the (percentage) rate fluctuations 
{Ax{t)} = {Ax} in G7 markets is described by an 
exponentially-decaying pdf p{Ax, At) where: Ax{t) — 
[xlt + At) - x{t)][100/x{t)]. This statistical behavior is 
not easily seen from the power spectrum of {x(t)}. Fig- 
ure ^ shows the p(Ax, Ai)'s for the Japanese 1^) and 
Philippine ^p) markets. For the G7 |h| and Philippine 
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FIG. 1: Asymmetric p( Ax, At) 's: a) Japan, and b) Philip- 
pines for At (day) = 1 (filled squares), 2 (diamonds), 5 (cir- 
cles), and 20 (cross-hairs). For any At, p(Ax,At) is nor- 
malized i.e., p{Ax, At)Sx — 1, where histogram bin-width 
Sx = (2/201)[max(Aa::niax, | — AxminI)]. Summation is for all 
Ax's. 



markets, our core data sets for {x{t)} comprise of 
daily closing rates from January 1971 to December 1993 
(M = 5000 trading days. At = 1 day). The Philippine 
market was under exchange control in the 70's and its 
p{Ax, At) has a large number of points around Ax — 0. 

Most p{Ax, At)'s exhibit a degree of asymmetry about 
Aa; = 0. An asymmetric p{Ax, At) is described by a pair 
of exponents ai and for Ax < (left) and Ax > 
(right), respectively. It is displayed by currencies that 
have either appreciated (e.g. Japan) or depreciated (e.g. 
Canada, Philippines) against the USD during the sam- 
pling period. 

Table HI lists the best-fit values (by least-squares 



2 



G X| 









Gain 






/° \l 


o 

e = 


N - 1 














Loss 









FIG. 2: Random walk in complex plane demarcated by an- 
gle O, into gain (+) and loss (— , shaded) regions. Each 
Aa; is the outcome of N random steps from origin O where: 
d = + a;f)^^^, and \Ax\ = cP . Sign of Aa; is given by 
location of endpoint N(a;fl, xi) in the plane. The n-th pha- 
sor is the fluctuation |Aa:„| = |a„p, of the n-th dealer where 
n = 1,2, ,iV. 



method) of a/ , and the left and right intercepts [3i and 
Pr of p(Aa;, Ai = 1 day) with the Aa; = hne for the G7 
markets. The exponential behavior of p( Ax, At) persists 
even for Af's > 2 days. We also found that p{Ax, At) 
broadens with At in a power-law trend. 

The behavior of {x{t)} has been described previously 
as a Markov process with multiplicative noise and the 
emergence of a highly-symmetric p{Ax, At) such as that 
of the German market 0, Q , is explained by the Fokker- 
Planck equation. However, these studies did not explain 
the emergence of asymmetric p{Ax, At)'s and the power- 
law behavior of their broadening with At. 

Here, we show that the statistics of {Aa;} is accurately 
described by a random walk on a demarcated complex 
plane. Our model accounts for possible asymmetry of 
p{Ax,At), its broadening with increasing At, and the 
power-law behavior of this broadening for G7 markets. 
Each Ax is the outcome of N random steps that start 
from the plane origin O (see Fig ^ . Its magnitude is 
Ax I — cP, where d is the Euclidean distance of the iV-th 
step endpoint N(xfl,,x/) from origin O, and Xfj and x/ 
are the real and imaginary coordinates. The sign of Ax 
depends on the location of N{xr, x/) in the plane which 
is demarcated by the polar angle O, into two regions 
representing the gain (-I-) and loss (— ) regimes. 

Random walk in complex plane. The location of 



N(xii:,X7) is given by the sum of N elementary pha- 
sors {an} H, N(xfl,x/) = Yln=i a«(a;_R,x/) = 
'^n=i kn|exp(j 9n) whcrc j = (—1)^/'^. The amplitude 
|a„ I and phase 0„ of the n-th phasor are statistically inde- 
pendent of each other and of the amplitudes and phases of 
all other {N — 1) a„'s. Possible |a„| values are uniformly 
distributed within a predetermined (non- negative) range. 
Phases 6n are also uniformly distributed: < 9n < 27r. 

Consistent with the central-limit theorem, the location 
of N(xii:, xi) is governed by a Gaussian pdf {N oo) 0: 
Pn{xr,xi) — (27rcr^)~^ exp[— (x|j -I- Xj)/2(t^], where ct^ 
is the variance. Hence, the rf^-values obey a negative- 
exponential statistics: Pd(d'^) — (2(7^)^^ exp(— (i^/2(T^). 
The q-th moment ((d^)?) of is given by: ((d^)') = 
q\{2a^)q = q\{d'^)'^, where the mean value (c?^) is 2a^. 
Phase 9 of N(xfl,x/) obeys a uniform statistics: peiO) = 
(27r)~^ for — TT < 9 < tt; pe{9) = 0, otherwise. 

The joint (second-order) pdf of d\ and d^ at 
two different time instants is j^: pd((i^,d2) — 
[{d^f{l - W)] exp {dl + dl) I [(d2)(i _ |^|2)] } 

Jo{2did2H/[(d')(l-H')]}, where: (d?) = (d^) = 
(d^), Jq{x) is the zero-order Bessel function of the 
first kind, and is the modulus of a complex factor 
that measures the degree of correlation between events 
at two different time-instants. If — (no correla- 
tion), pd(df,d2) — Pd{dl) pdCdj). On the other hand, 
Pdidl,dl) = Pd{d\) 5{d\ - dl), as \fi\ 1, where S{x) is 
the Dirac delta function. 

As it is, the statistics of a random walk in a complex 
plane is insufficient to describe the characteristics of the 
p(Ax, At)'s of forex markets because the possible Ax's 
are from — oo to oo, while d^ > 0. 

Demarcated complex plane. We solve the problem by 
demarcating the complex plane into two sectors repre- 
senting the gain (-|-Ax) and loss (—Ax) regimes, where 
we identify that: |Ax| = d^. The gain area is set by the 
polar angle O (counterclockwise rotation) and Ax is pos- 
itive (negative) if N(x/j,x/) is on the gain (loss) sector. 
The gain (loss) area is zero if O = (O = 27r). 

Figure|31presents exponentially-decaying histograms of 
Ax's generated at 8 = 0, 7r/2, tt, 37r/2, and 27r (360 
deg). The histograms are asymmetric about Ax = at 
all angles except at O = tt where the plane is equally 
divided between the gain and loss sectors. The corre- 
sponding 8's for the G7 markets in Tableware not easily 
determined since Q is not clearly related with ai, ar, Pi, 
and Pr- However, our model reveals a unique relation 
between R = Ai/Ar and 8, where Ai (Ar) is the area 
under the best-fit exponential for Ax < (Ax > 0) of 
the histogram. 

Figures 0^-b plot the average (R) as a function N and 
number of trials, respectively (8 = tt). (R) is insensitive 
to N but for a fixed N, the standard deviation of (i?) 
decreases quickly (power-law decay) with the number of 
trials . Figure 21: shows that the dependence of (R) with 
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TABLE I: Parameter values of exponential fit: /3 exp(aAa::) 
p{Ax,At): ai = \ar\, and ln{/3i) =ln(|/3r|). 



to p{Ax, 1 day) for G7 and Philippine markets. For a symmetric 



Country 



ln(A 



Canada 
Prance 
Germany 
Italy 
Japan 
UK 
Philippines 



6.143 ±0.197 
2.283 ± 0.082 
2.166 ±0.076 
2.316 ±0.113 
2.371 ±0.115 
2.171 ±0.091 
4.711 ±0.904 



-5.807 ±0.150 
-2.469 ± 0.096 
-2.327 ± 0.068 
-2.44 ±0.102 
-2.635 ± 0.084 
-2.244 ± 0.086 
-2.892 ± 1.063 



1.090 ±0.071 
0.056 ± 0.065 
0.058 ± 0.062 

-0.012 ± 0.090 
0.073 ± 0.092 
0.004 ± 0.074 

-1.745 ±0.564 



1.037 ±0.058 
0.181 ±0.069 
0.125 ±0.050 
0.152 ± 0.074 
0.187 ±0.062 
0.063 ± 0.066 
-1.464 ±0.964 
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FIG. 3: Random walk in demarcated complex plane. His- 
tograms of Aa;s for = (squares), n/2 (diamonds), tt 
(circles), 37r/2 (triangles) and 2n (crosses) where N — 50, 
< [a„| < 0.015 (double precision). 



e is well-described by: (i?) = 6/(360 dcg - 6). 

The (i?)'s are calculated from p(Aa;, Ai)'s of the G7 
markets and then used to determine their correspond- 
ing 0's via the i?(9)-curve. The following 8's (deg) 
were obtained (At — 1 day): Canada (179.7), France 
(180.4), Germany (182.3), Italy (177.8), Japan (177.6), 
and UK (174.9). For the G7 markets, differences in the 
size of the gain and loss sectors are small (< 4%). On 
the other hand, the p{A.x, 1 day) for the Philippine mar- 
ket is highly asymmetric with a significantly small gain 
region (0 = 108.1 deg). Against the USD, the Philippine 
peso depreciated by 246% from December 1983 (lifting of 
exchange control) to December 1993. 

Broadening of p{A.x, At) with time delay At. For 
the G7 markets, p{Ax, At) broadens with increasing At, 
while preserving its original negative exponential statis- 
tics within: 1 < At (day) < 20. Previous studies on 
broadening were confined to At < 12 hours 

FigurcElplots the dependence of ai and ar with At for 
the G7 markets where a power-law behavior is observed. 
Except for Canada, the dependence of aj's (and jail's) is 
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FIG. 4: {R) as a function of (9 = tt): a) N, b) Number of 
trials per datapoint (A'^ = 50) ((s.d.) = 1.05.), and c) Q for 
N — 50. 10 trials for each datapoint. 



remarkably described by one and the same power-decay 
curve (Reference curve 2) indicating a scale-free behav- 
ior for the market dynamics. For the Philippine market, 
the At-dependence of the ai (ctr )-values is erratic with 
large standard deviations - there is increasing asymme- 
try for p{Ax, At) with At. Exchange control destroyed 
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FIG. 5: Dependence of ai (positive) and Or (negative) 
with At for Canada (squares for oj, solid squares for Qr), 
France (right triangles), Germany (inverted triangles), Italy 
(diamonds), Japan (triangles), UK (circles) and Philippines 
(crosses for crossed squares for Ur)- Reference curves: 
ai = 6.4At- 0.561 (Curve 1) and a; = 2.21At- 0.593 (Curve 
2). p(AxQ, Q) contains 5000 datapoints and each point is the 
average of 10 trials. 



the scale-free behavior of the Philippine market. The 
largest standard deviation for Canada and UK is ±0.2 
and ±0.1, respectively which happens at Ai — 1 day, 
while for Philippines it is ±1.4 (At — 3 days). For 
any G7 market, (Az) increases with At according to a 
power law since the exponent is inversely proportional to 
1/2^2 = l/(Ax). 

The observed power-law dependence of the broaden- 
ing of p(Aa;, At) with At is explained as follows. Let 
(Axq) = {dr)^, be the average fluctuation correspond- 
ing to a longer time-delay Atg = QAt, where: dr = 
\E,d,eMjO,)\ = [(E,rf,cos0,)2 + (E,d,sin0,)2]i/2^ 
At is the basic time delay, and q — 1,2, . . . ,Q. If (di) = 
{d^) = ... = (dg), then (dr) = (di) ([(E, cos 0,)^ + 
(E,sin0,)2]V2)^QP(rf^). 

The (dx) value depends on the probability distribution 
of {Oq}. If the 9q^s are statistically independent and uni- 
formly distributed within [0, 27r] then dx ~ Q'^^^{di}, i.e. 
p « 0.5. Figure ini inset plots ai and la^l of p{Axq,Q) 
with Q for the above case which verifies our assumption 
of a power-law dependence of {dr) with (di). 

Figure ini plots p as a function of uncertainty spread SO. 
If 0g's are restricted such that {q > 1) is uniformly- 
random within the forward (— 7r/2 < 9q — 0.566 < Oq+i < 
Oq + 0.566 < 7r/2) and opposite (7r/2 < {Oq + Tr) -0.569 < 
Oq+i < {Oq + tt) + 0.556' < 37r/2) directions then: < 
p < 0.5. If the Oq^s occur only within the forward or 
opposite direction then 0.5 < p < 1. In all cases, 9i is 
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FIG. 6: Power p vs. Possible directions of O^'s are always 
forward or always backward (filled circles). Each datapoint 
is average of 10 trials. Reference curve: p — (iff)'^''^'' . Inset: 
Of; and a,, of ln[p(Aa;Q, Q)] vs Q for Q (deg) = 90 (squares), 
180 (circles), and 270 (diamonds). Reference curve: a; = 
76.6Q~\ 



uniformly-random in the range [0, 27r]. 

A G7 market exhibits a (At)~°-^-dependence (i.e. p w 
1/4) of Oil (Fig|SJ) because dg+i is equally likely to be in 
the same or opposite direction of dq, with 6Q « 0.3 (17 
deg). For a market where the directions of dg+i and dq 
are statistically independent of each other, the decay of 
a/ with At is faster (p « 1/2, FigElinset). 

In our model, the forex market consists of N indepen- 
dent players where Ax is treated as a random walk of N 
steps in the demarcated complex plane where |a„p rep- 
resents the fluctuation |Aa;„| of the 7i-th player. Plane 
anisotropy leads to asymmetric prf(c?^)'s. Our model 
also explains the power-law dependent broadening of the 
p{Ax, At)'s with At in G7 markets. 

Ax is interpreted as the intensity of the resultant com- 
plex amplitude that arises from the linear superposition 
of the fluctuations of the N independent dealers. In- 
teractions between players arises inherently because the 
intensity which is the product of the resultant complex 
amplitude and its conjugate, contains interference terms 
between the contributions of the individual players. We 
showed that real market dynamics could be analyzed ac- 
curately with a relatively low number of the interacting 
players. The interaction between the multitude of agents 
in a real market could be effectively reduced into one with 
a low number of representative classes. 
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Our model currently neglects the phenomenon of al- 
lelomimesis which causes social agents to cluster [loj. 
Herding (bandwagoning) among market dealers could 
trigger massive selling or buying in financial markets and 
leads to large swings in {x{t)} like those found in the 
middle 1980's for the G7 markets. The same limitation 
is found in previous models of forex markets 0, Q . 
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